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We study a supremacy distribution in evolving Barabasi-Albert networks. The supremacy Si of 
a node i is defined as a total number of all nodes that are younger than i and can be connected to 
it by a directed path. For a network with a characteristic parameter m = 1, 2, 3, . . . the supremacy 
of an individual node increases with the network age as in an appropriate scaling region. 

It follows that there is a relation s{k) ~ fc™^^ between a node degree k and its supremacy s 
and the supremacy distribution P{s) scales as s~i~^/(i+'"' . Analytic calculations basing on a 
continuum theory of supremacy evolution and on a corresponding rate equation have been confirmed 
by numerical simulations. 



PACS numbers: 



i.75.-k, 02.50.-r 



I. INTRODUCTION 

During the last few years there has been a large in- 
terest in modeling of networks 0, 0, 0, S 13 and several 
parameters describing the network structure have been 
considered. The examples arc: degree distribution P{k) 
[2, iS] , mean path length (7,, ^ ^ J^] , betweenness cen- 
trality fload) M . [T ill or first and higher order clustering 
coefficients Il3l Il4) . Universal scaling has been ob- 
served for some of these parameters in computer simula- 
tions and in real data describing such objects as the In- 
ternet, WWW, scientific collaboration networks or food 
webs 1^ 0, 01 . Here we introduce a new parameter that 
can play an important role for description of a class of 
directed networks. We name the parameter a supremacy 
since it describes the number of nodes that are subordi- 
nated to a certain node. In the next Section we define 
our parameter and show its relevance for different prob- 
lems of complex networks, Sec. IIIII includes a continuum 
theory for the supremacy time evolution Si{t) and the 
supremacy probability distribution P{s) in the Barabasi- 
Albert (BA) model with m = 1, in Sec. IIVI we find and 
solve a corresponding rate equation while in Sec. 0a gen- 
eralization of our problem for the BA model with m > 1 
is presented. 



II. THE MODEL 

Let us consider the BA network with the characteris- 
tic parameter m = 1 P, 0] . At the moment ti a node 
i is created and it attaches to some older node in the 
network according to the preferential attachment rule 
(PAR). Then in the next time steps other nodes are cre- 
ated and are attached to the node i or to other nodes of 
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FIG. 1: Schematic illustration of the supremacy effects in 
tree- like BA network with m — 1. Numbers situated in the 
vicinity of the nodes represent their supremacies. 



the network following PAR. As a result at the moment 
t > ti there is a subgraph in a form of tree T{i,t) be- 
ginning in the node i and containing all nodes that are 
younger than the node i and that are connected to i by 
directed paths as in Fig. ^ If we assume that the node 
i represents a scientist who wrote an important paper or 
a politician who created an influential party [l^ UM 
can consider all nodes belonging to the tree as his/her 
successors. If the tree T{i,t) contains Si nodes then the 
number is the measure of the supremacy or the pre- 
dominance of the node i at time t. The subgraph T(i) 
can be also interpreted as a cluster of connected sites in 
the directed percolation problem [Tl[l9ll20l l2lll2^ 
and the supremacy of a node i is just the size of such 
a cluster starting from the site i. Since the evolution of 
the network is governed by PAR and all properties of the 
network are described by some probability distributions 
we are interested in the supremacy distribution P{s) in 
the network. 
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FIG. 2: Supremacy as a function the node degree. Solid lines FIG. 3: Supremacy distribution in BA model. Solid lines 
represent analytical predictions of s{k) given by Eq|Klin the represent analytical predictions of P{s) given by EqElin the 
case of m = 1 and Eq^jin the case of m = 2. case of m = 1 and Eq^]in the case of m = 3, 5. 



III. CONTINUUM THEORY OF SUPREMACY 
EVOLUTION AND DISTRIBUTION FOR m = 1 

To find the supremacy distribution P{s) we follow the 
method that was introduced in 2] for calculation of de- 
gree distribution P{k) in evolving networks. We start 
from determining the time dependence of Si (t) assuming 
that it is a continuous real variable. The supremacy of 
the node i increases in time because new nodes can be at- 
tached to any node of the tree T{i,t). Let nodes belong- 



ing to the tree T{i,t) possess degrees k^j^\kf'\ 



Using the PAR we can write the following equation for 
changes of Si{t) 



ds,{t) 
dt 



(0 



^ 2t 



1=1 



2t ' 



(1) 



where K{i) — J2tLi ^-i^^ used the fact that at 
the moment t the sum of all nodes degrees in the whole 
network equals to 2t. On the other hand taking into 
account the tree structure of the considered subgraph we 
can write the supremacy Si as 



1=1 



thus K{i) — 2si 



with the solution 



1 and we have a simple equation 

ds,{t) _2s,~l 



dt 



2t 



S^{t) 



1 



(2) 



(3) 



(4) 



where we took into account the initial condition Si{t = 
ti) = 1. The solution Q means that the node supremacy 



increases linearly in time comparing to the square root 
dependence of the node degree 0, i.e. ki{t) = ^J~^- 
Combining the last two results we get a simple relation 
between the node supremacy and the node degree 



m = \[k^ + 1] 



(5) 



In the region k < 100 this formula fits well to numeri- 
cal simulations presented in FigI21 while for larger k dif- 
ferences between the analytic theory and the numerical 
simulations are observed. 

The probability density P{s) for the supremacy distri- 
bution in the network follows from the relation 



(6) 



where P{ti) = l/< is the distribution of nodes attachment 
times ti for a network of age t. After a simple algebra we 
get 



{2s, - 1) 



(7) 



One can see that the supremacy distribution is a time 
independent function. Fig|31shows the comparison of the 
last equation to numerical data. Let us stress that for 
s » 1 the supremacy distribution scales as P{s] ~ 
while the degree distribution for BA model P, Q scales 
as P[k) - /fc-3. 



IV. RATE-EQUATION FOR SUPREMACY 
DISTRIBUTION FOR m = 1 



Now we show how to get the supremacy distribution 
using the rate-equation approach that was introduced by 
Krapivsky, Redner and Leyvraz 6] to study networks 
degree distribution P{k). Let N{s,t) is the number of 
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the supremacy s at time t. The rate 



nodes possessing 
equation for N{s,t) is 

dN{s, t) _ [2{s - 1) - l]N{s - 1, t) - (2s - l)N{s, t) 



dt 



2t 



+Ss,i.{S) 



The first term on the right-hand side of Eq. © corre- 
sponds to creation of a new node with the supremacy 
s. The process is proportional to the number of nodes 
with the supremacy s — 1 and the corresponding transi- 
tion probability that follows from the PAR and Eq. ij^J. 
The second term corresponds to creation of a node with 
a supremacy s + 1, i.e. to destruction of a node with a 
supremacy s while the last term describes creation of a 
node with a supremacy s = 1. Writing N{s, t) — P{s)Nq 
where TVq = t corresponds to the total number of nodes 
at time t and P{s) is the probability of a node with the 
supremacy value s we get the recursive equation 



P{s) 



2s 



-P{s - 1) for s > 2, 



2s + 1 

where P(l) = 2/3. The solution of Eq. ® is 
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Pis) 



(2s-l)(2s + l)' 



(9) 



(10) 



Note that for s ^ 1 the solution H10|l coincides with the 
solution that has been received in the limit of the 
continuum theory. 



V. SCALING OF SUPREMACY 
DISTRIBUTION FOR m > 1 

The peculiar feature of the BA model is the indepen- 
dence of the scaling exponent characterizing the degree 
distribution P(fc) ~ from the model parameter m 
describing the number of links that are created by every 
new node. Below we show that the scaling exponent of 
supremacy distribution depends on the parameter to. If 
we neglect all loops existing in the BA network with the 
characteristic parameter m > 1 then we can easy repeat 
our considerations from Sec lIIII and IIVI Instead of Eq. 
@ we get 



s^ = H-^(fcf -TO) = is:(i)-M-TOS„ (11) 
1=1 

and time evolution of the supremacy is described by 



TO it 



TO + 1 V ti 



1 



TO + 1 ' 



(12) 



thus the relation between the degree and the supremacy 
is 



, , \ m+l -I 

m k \ 1 



TO + 1 V TO 



TO -|- 1 



(13) 




FIG. 4: Schematic illustration of supremacy effects in BA net- 
work with m = 2. Solid arrows represent connections within 
the supremacy area / cluster of the top vertex, whereas dashed 
arrows express connections pointing outside the cluster. Note 
that there is a single loop in the cluster. 



It follows that for dense networks with to 3> 1 the 
supremacy Si{t) increases in time much faster than the 
node degree ki{t). Fig. [3 shows a comparison of the re- 
sult (|13|l to numerical data for to = 2. One can see that 
the predicted scaling of s(fc) breaks down completely for 
large values of k where the plot s(fc) saturate. The rea- 
son is the presence of loops that for to > 1 appear in the 
network and that have been neglected in our approach. 
If TO > 1 the result \i'A\ is valid mainly for vertices with 
a small degree ki (and a small supremacy Si) since loops 
are sparse in small clusters starting from such nodes. The 
saturation effect does not appear for the BA model with 
the parameter to = 1 where loops are absent. 

Taking into account (|12|l we get the supremacy distri- 
bution in the form 



P{S^ 



1 


dsi 


' 2 


(to + l)Si — 1 


1 


dt. 


TO 


TO 



^ + 3 
1 + 1 



(14) 



We see that the scaling exponent for the supremacy dis- 
tribution equals to5 = —l — 2/(1-1- to) and in contrast to 
the scaling exponent of degree distribution it depends on 
the parameter to. The result (|14|l is in a good agreement 
with numerical simulation for BA networks, see Fig. |3| 
The rate-equation for to > 1 is similar to Eq. ((SJ, i.e. 

dN{s,t) [{l + m){s-l) -l]N{s-l,t) 



dt 2t 

[il + m)s-l)]N{s,t) 
2t 



+ Ss,i (15) 



The resulting solution for the probability P{s) can be 
written as the following product 



Pis) 



TO 



n 



[iim+1) + 1] 



(16) 
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for s > 1 where P(l) = 2/(m + 2). For dense networks 
TO 3> 1 the sohition can be approximately written 
as 



Pis) 



ms 



(17) 



what coincides with (|14|l . 



VI. CONCLUSIONS 



follows that there are universal scaling laws describing 
the time evolution of the supremacy and corresponding 
supremacy distributions in BA models. On the contrary 
to the scaling results for nodes degree the corresponding 
scaling exponents of supremacy depend on the charac- 
teristic model parameter m. Numerical simulations are 
in good agreement with analytical estimations for node 
with a small and medium supremacy especially for the 
case TO = 1 where no loops are present in the system. 



In conclusion, we introduced a universal parameter (a 
supremacy) that describes vertices in directed networks. 
The parameter equals to the size of a cluster starting 
from the site in a directed percolation model. We have 
shown that for Barabasi- Albert model there is a rela- 
tion between the supremacy and the vertex degree. It 
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